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We argue by observation of the neural data that neurons in area dMEC of rats, which fire whenever the 
rat is on any vertex of a regular triangular lattice that tiles 2-d space, may be using an advanced numeral 
system to reversibly encode rat position. We interpret measured dMEC properties within the framework of 
a residue number system (RNS), and describe how RNS encoding - which breaks the non-periodic variable 
of rat position into a set of narrowly distributed periodic variables - allows a small set of cells to compactly 
represent and efficiently update rat position with high resolution over a large range. We show that the uniquely 
useful properties of RNS encoding still hold when the encoded and encoding quantities are relaxed to be real 
numbers with built-in uncertainties, and provide a numerical and functional estimate of the range and resolution 
of rat positions that can be uniquely encoded in dMEC. The use of a compact, 'arithmetic-friendly' numeral 
system to encode a metric variable, as we propose is happening in dMEC, is qualitatively different from all 
previously identified examples of coding in the brain. We discuss the numerous neurobiological implications 
and predictions of our hypothesis. 



Recent experiments reveal that rats, notoriously good 
navigators, encode information about their position in 2- 
dimensional spaces in a remarkable way Ol|^. Each neuron 
in the dorsolateral band of the mediolateral entorhinal cortex 
(dMEC) fires when the rat is on any vertex of an imagined 
regular triangular lattice, tiling the plane. The firing pattern 
is independent of the enclosure size and shape, and updates 
correctly as the rat moves around, even in complete dark- 
ness. These observations, together with lesion studies, hint 
that dMEC may play a central role in rodent path integration 

SiBi. 

Nearby dMEC neurons share the same lattice period (k) and 
orientation, and differ only in their relative spatial phases 0]. 
Thus at any instant, the active subset of neurons only specifies 
cuiTent rat position as a phase: upto, or modulo, all possible 
periodic displacements of the lattice. Moving ventrally along 
the length of dMEC, the neural lattice period increases mono- 
tonically, but the measured range is naiTow: from approxi- 
mately 30 cm to 70 cm J3l. In summary, we note that dMEC 
neurons decompose the 2-dimensional vector of rat position 
into a set of phases, modulo a set of lattice periods, along two 
independent lattice directions. 

While phases within the smallest lattice - about the size 
of the rat - might arguably vary at the appropriate scale for 
fine position discrimination during navigation, the largest lat- 
tice period (<2 m, by extrapolation along dMEC length) puz- 
zlingly appears to fall far short of the range over which a rat 
might know its position. 

Can the phases of such a narrow range of lattices even theo- 
retically contain enough information to unambiguously repre- 
sent position over the behaviorally relevant range of distances 
covered by rats? And even if so, what is the advantage of 
decomposing rat position - a quantity that can in theory be 
represented simply by a pair of numbers - in this distributed, 
seemingly bizarre way? 

Based on the observation that dMEC stores only mod- 



ulo information about rat position, we suggest that dMEC 
may be encoding and enabling the reconstruction of the two- 
coordinate rat position vector according to a generalized ver- 
sion of a residue number system (RNS) |6], a scheme with 
uniquely useful properties in the neurobiological context. In 
an RNS, a number x (e.g., position in 1-d) is represented by 
a list {xi = X mod Xf) of its residues, or remainders after divi- 
sion by a set of fixed numbers (A-i , A-2 , . . . , A-at), relatively prime 
to each other, called moduli. For example, if the moduli are 
(13, 15, 16, 17, 19), the number 1000000 is represented by the 
residues (1, 10, 0, 9, 11). The Chinese Remainder Theorem 
(CRT) guarantees that any number smaller than the product of 
all the moduli is uniquely specified by (and therefore can be 
reconstructed from) its residues. In the example above, any 
number up to 1007760 has a unique representation. An RNS 
has striking parallels to dMEC encoding: According to our 
interpretation, the dMEC lattice periods are the RNS moduli 
(which need not be co-prime integers, as we discuss later); 
the lattice phases, specified by the active set of neurons at any 
given rat position, are the residues. 

In an RNS, large numbers (e.g. position in a large space) are 
represented combinatorially and thus compactly; this prop- 
erty is shared with base numeral systems like decimal or bi- 
nary. A small set of registers (A^) can represent a very large 
range of numbers ([0, 10^ — 1]: decimal, [0,2^ — 1]: binary, 
or [0, A,^ — 1]: RNS, where X is the approximate size of 
the moduli, and is the number of moduli; in the example 
above. A, « 16, = 5). By contrast, in a sparse encoding 
scheme where each element represents one possible number 
(e.g. place-cell representation of position in hippocampus), 

elements only cover the range [0,A^— 1]. But besides its 
combinatorial capacity, two unique features of an RNS make 
it especially useful, compared to base numeral and other sys- 
tems, for position encoding in dMEC. 

First, the set of registers in an RNS that support its vast 
representational capacity can be very closely spaced, span- 
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ning a narrow range: as illustrated above, all moduli may be 
nearly equal. By contrast, successive registers in base numeral 
systems must represent information on geometrically spaced 
scales differing by powers of say 2 or 10. The largest register 
must represent a scale comparable to the largest representable 
number and the smallest register scale must be comparable 
to the finest resolvable detail. The ability of an RNS to use 
closely spaced registers to represent a large range of positions 
is critically useful for position encoding in dMEC, with its 
narrowly distributed lattice periods. 

Second, addition, subtraction, and multiplication are com- 
pletely parallelized in an RNS, because unlike base numeral 
systems, they require no "carrying-over" of information from 
one register to the next: the sum of two numbers is the mod- 
ulo sum of their residues computed independently within each 
register In computer science, this property has long been ap- 
preciated and is used to perform fast, parallelized computa- 
tions in signal processing and RSA applications The 
moving rat, using an RNS, can update its estimate of posi- 
tion by independently incrementing the phases of each lattice, 
without carrying phase winding information from one lattice 
to another Indeed, in dMEC each lattice phase is indepen- 
dently updated as a function of rat position, without evidence 
of jumps in phase when a smaller lattice completes one wind- 
ing. 

Some important properties of dMEC encoding differ from 
a standard RNS. Lattice periods in dMEC are probably not 
co-prime integers, but real numbers, as is rat position. Also 
unlike an RNS, where the number of residues (distinguishable 
phases or phase resolution) grows with modulus size (lattice 
period), the phase resolution in dMEC seems constant across 
lattices. Under these conditions, existing exact formulas for 
reconstructing numbers from their residues using the CRT no 
longer apply, and the representational capacity may be much 
smaller than possible in a standard RNS. But despite these 
differences, summation and multiplication are still fully par- 
allelized, and we show without recourse to any particular de- 
coding scheme, that the representational capacity continues 
to scale exponentially with the number of (non co-prime) lat- 
tices, spanning a vast range of real-valued positions with high 
resolution (Figure 1 and Supplementary Information). The 
uniqueness of this representation guarantees that an inverse 
map from lattice phases to position still exists. The capac- 
ity (Figure 1) under a generalized RNS encoding is far in ex- 
cess of the navigational requirements of rats, and the surplus 
capacity could be devoted to redundancy for error correction 
^,8J and robustness. 

To summarize our argument, we suggest that dMEC treats 
position simply as a 2-coordinate vector, and then encodes the 
coordinates, two real numbers, using a unique numeral system 
that can represent large numbers with a narrowly spaced set of 
registers. At first glance, one might have imagined that dMEC 
firing reflects some kind of integral transform (e.g. Fourier or 
Wavelet) of 2-dimensional space and the rat's location in it. 
To do this, dMEC would have to treat rat position as an entire 
2-d function, such as a blob centered at the rat's location in 
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FIG. 1 : Capacity of dMEC under generalized RNS encoding. The 

maximum uniquely representable distance, D, grows exponentially 
as a function of number of lattices, Nx, a, and as a power of the 
phase uncertainty A(|), b, even when the lattice periods are not co- 
prime and even if the phase uncertainty is equal across lattices. [The 
first lattice period is 30 cm, with 4 cm increments per subsequent 
lattice. In a A(|)/27r = 0.2, and in b A^x = 12. The fit parameters are 
a ~ 0.55 and ^gff— 9.7. See Supplementary Information for more 
details.] c 12 lattices uniformly spaced from 30 cm to 74 cm with 
5 resolvable phases in each direction can unambiguously represent a 
•-^ 2 km X 2 km area with 6 cm resolution. If 5000 neurons build 
each lattice, that would require ~ 5 x 10'* neurons. To cover the 
area with sparse unimodal place-cell like encoding (with 10 neurons 
per (6 cm)^ block) would require ~ lO"^ neurons, compared to the 
estimated 10^ neurons in rat dMEC 1 16, 17], which could sparsely 
represent at most ~ 6 mx6 m with (6 cm)^ resolution. With 24 
lattices, £) ~ 2 x 10^ km in each direction, hugely in excess of the 
representational requirements of rats. 



space, and transform this into a set of 2-d patterns of different 
spatial frequency. Besides the computational unnecessity of 
representing two position coordinates by an entire function, 
the integral transform scheme is unlikely for at least two rea- 
sons: (1) It is incompatible with the data. The experiments 
iB Q unequivocally show that the patterns of dMEC ac- 
tivity are stereotyped, and independent of enclosure size or 
shape; given the activity pattern of one neuron in each lattice 
in any one enclosure, the activities of the rest are well speci- 
fied upto a spatial phase shift in all enclosures. Thus, dMEC 
encoding contains little additional information or flexibility to 
describe general functions in 2-d space, as required of inte- 
gral transforms. (2) A Fourier-Uke transform cannot be used 
to represent position unambiguously over a distance greater 
than the largest lattice period (^ 2 m). To reconstruct a sin- 
gle bump of rat position over a range D with resolution Dmin 
would require lattice periods ranging from smaller than D^m 
up to D. 

Our focus here has been on characterizing the general the- 
oretical properties of position encoding in dMEC. These in- 
sights provide the necessary foundation for understanding 
how information encoded in dMEC may be used by areas 
that receive dMEC inputs. We have illustrated, without re- 
gard to specific decoding schemes, how the maximal capacity 
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of dMEC for unambiguous position encoding scales as a func- 
tion of dMEC phase resolution and lattice numbers, and have 
observed based on the properties of an RNS and our simula- 
tions, that these phases theoretically contain enough informa- 
tion to uniquely encode position over a combinatorially large 
range. As we discuss next, our interpretation of dMEC encod- 
ing raises several testable questions with specific implications 
for theory and experiment. 

The most basic of these is whether the rat does actually 
make use of the range of information we have shown is theo- 
retically available in dMEC under an RNS -like scheme. There 
are two (not mutually exclusive) ways in which rats may use 
this information: (1) For homing and path integration over 
large ranges by decoding the relative dMEC phases to com- 
pute rat displacements, or (2) For attaching unique 'labels', as 
described below, to a large number of specific locations, with 
landmark-independent path integration only between nearby 
locations. Either scenario can be probed experimentally, as 
we describe below, but in both cases resolving the fundamen- 
tal issue requires testing in enclosures of size approaching the 
behaviorally relevant range for rats. 

In the first case, assuming dMEC is the primary source of 
position representation, our proposal would be strongly sup- 
ported if the rat can perform reasonably accurate landmark- 
free homing behaviors over ranges much larger than the 
largest dMEC lattice period. This scenario would require 
an explicit decoding of the phase encoding of position. Be- 
cause position encoding by modulo residues, which we argue 
is happening in dMEC, is a 1-1 onto (bijective) function over 
the illustrated range (Figure 1) of positions, the mathemat- 
ical inverse exists and is unique. However, there could be 
numerous possible neurobiological decoding schemes to ex- 
actly implement this inverse or approximate it. Our capacity 
estimate, which contains the experimentally observed phase 
uncertainty, bounds how much information can be extracted 
from dMEC and should be used as a guide for evaluating the 
efficiency of different decoding schemes. Experimentally, the 
test of the long-range position encoding hypothesis involves 
measuring the largest lattice period in dMEC, and quantita- 
tively determining the range over which rats can perform ac- 
curate homing in large featureless spaces. 

In the second case, the large set of unique dMEC phases 
as a function of rat position may be used as absolute markers 
for specific landmarks or positions in familiar environments. 
In this scenario, the rat may not explicitly decode position or 
displacements over distances larger than a lattice period, but 
may essentially use the vast set of distinct dMEC phases to 
uniquely represent a large number of distinct, spatially sepa- 
rated locations. The rat may locally perform path integration 
by updating phases, to ensure that starting from an absolute 
phase at landmark A, the rat can take any path to landmark B 
and yet obtain the correct absolute phase for B. Such a system 
would be useful to the rat in distinguishing a familiar land- 
mark from a lookalike but spatially distinct location fioll . This 
proposal could be tested experimentally by checking at multi- 
ple familiar landmarks in large, partially occluded enclosures 



whether the absolute dMEC phases for each landmark are re- 
producible across trials. If dMEC is reset to a set of previously 
assigned phases for landmark B upon reaching B from a novel 
occluded path and after extensive exploration in an unfamil- 
iar landmark-rich environment, it would suggest that absolute 
phases are important; the RNS scheme would allow dMEC to 
uniquely encode a large number of landmarks. 

Unlike most numeral systems, in an RNS scheme the repre- 
sentation of space by the different registers or dMEC lattices 
is not hierarchically ordered: all lattices are roughly equal 
in their contribution to position representation at any scale. 
Therefore, if hippocampal place cells are involved in recon- 
structing a position estimate or reading out a unique label of a 
landmark from dMEC, each local group of hippocampal cells 
must pool input from many different dMEC lattices; indeed, 
neurons in the septal (dorsal) half of the long axis of the hip- 
pocampus do appear to receive inputs from across the dMEC 
band 1 1 1 ]. In addition, position reconstruction from an RNS 
is susceptible to characteristic deficits following microlesions 
that successively destroy individual dMEC lattices. Error 
correction and redundancy could grant the system resilience 
against partial lesions, but beyond a critical point, there should 
be a sudden degradation in position estimation or location rep- 
resentation at all scales. By contrast, in any scheme where 
different lattices encode positions at their corresponding spa- 
tial scales, selectively lesioning the smallest (largest) lattices 
is likely to destroy fine (coarse), but not coarse (fine), position 
estimation. 

When rats are trained to explore small enclosures or 1-d 
tracks, neurons in the hippocampus - the primary output target 
of dMEC, and a critical locus of spatial learning and memory 
- tend to form unimodal place-fields that resolve position with 
relatively high accuracy, and cover the space ImlU. Based 
on our analysis, we predict that the inefficiency of sparse en- 
coding (Figure Ic) and limited neuron numbers in hippocam- 
pus (/^ 10^ — 10^ neurons) fliil . rule out the possibility that 
hippocampus could fully remap large enclosures with high 
resolution using unimodal place-cells. Rather, driven by rep- 
resentational limitations in large areas, place cells must either 
generically develop multi-peaked responses as a function of 
rat position to cover space with reasonably high resolution, 
akin to dMEC neurons, or must disengage from finely repre- 
senting all space with narrow place fields, instead covering 
only select locations based on other salient cues or associ- 
ations. The latter scenario would make hippocampus a se- 
lective consumer and processor of position information from 
dMEC, leaving it free to perform more general associative 
tasks ii[li[ll. 

The possibility that dMEC may be representing position, a 
continuous metric variable, using a compact and parallelized 
numeral system amenable to arithmetic operations such as ad- 
dition or shifts on the variable, is itself extraordinary. Such a 
numeral system code is qualitatively different from all other 
known examples of coding in the brain |18]: Proportional- 
rate coding (e.g., eye position and firing rate in the oculomo- 
tor system 1 19]), unary or sparse coding (e.g., head direction 
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cells fid], place cells fl2il '). or even other kinds of combina- 
torial codes which represent non-metric variables (e.g. odor 
iilll ') and therefore lack an arithmetic aspect. This encoding 
scheme provides insight into the ingenuity of neural codes, 
and provokes questions at the intersection of neuroscience, 
mathematics, and computer science. 
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We numerically estimate the range of positions that the 
phases of a small set of lattices, with spacings distributed 
over a narrow range, can unambiguously represent. We 
assume that unlike a standard RNS, the dMEC lattice 
spacings are not co-prime integers, and the encoded po- 
sition is a real number. There is no clear analytical state- 
ment for how the representational capacity of a general- 
ized RNS, with limited, constant phase resolution across 
lattices and non co-prime lattice spacings, compares to 
that of a standard RNS, but we show without regard to 
any particular readout scheme that the capacity is large, 
and scales similarly to a standard RNS. For simplicity 
and clarity of exposition, most computations are in one 
dimension, but can be extended to two dimensions. 

Exact numerical evaluation Consider a set of lat- 
tice spacings or modulii, — 1,...,N\}. A position 
xi is represented by the vector of its residues, 

xii = a;i mod Ai, i^l,...,Nx, 

or phase vector (j)i , whose entries are 



distinct states of the N\ registers. 



A. 



mod 2n. 



For example, when xi xi + A^, then (pi}^ {(pik + 
27r) mod 27r — (jji^.. Assume that each phase can only be 
resolved upto a spread A(j). Then the lattice phases can 
be used to distinguish between xi and X2 if there is a 
mismatch of more than A(j) in at least one entry of their 
phase vectors. 

Therefore, we define the phase distance between xi 
and X2 to be the largest phase mismatch across lattices 
(entries) in their phase vectors: 

II f^i - '?5'2 11= maxi[min[ {(f>u - (f>2t) mod 2Tr, (1) 

(27r - -I- ) mod 2tt]] 

In the numerics for our figure, the smallest lattice spac- 
ing is 30 cm, and each subsequent lattice spacing is 4 cm 
greater than the last. We conservatively assume a phase 
resolution of 1/5 of a period, or A(p — 2tt/5. We start 
with (pi = for a; = 0, and increase x in small incre- 
ments to generate Figures la and lb. The range D in 
Figure la is defined as the point at which the phase dis- 
tance between a; = and x = D drops to within the 
phase resolution, Acj). Figure lb plots the phase resolu- 
tion needed (or maximum allowed phase uncertainty) to 
cover the range a; = Otoa; = Z?asa function of D, for a 
fixed number of lattices {N\ = 12). 

Scaling ansatz for generalized RNS In a standard 
RNS, the range of sequential numbers that can be dis- 
tinguished from each other is bounded by the number of 
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The representable range reaches this upper bound if 
the registers are co-prime integers (Otherwise the rep- 
resentable range is the least common multiple of the reg- 
isters). To speak of distances instead of integer numbers, 
we may associate a length scale of unity to the finest dis- 
tance resolution, -Dmin = 1- The maximum representable 
range is then Dmax = A^max, as given by Eq. JSJ. If all 
the lattice periods A^ are similar in size, then the maxi- 
mum range scales exponentially with Nx and as a power 
of the distinct number of states in each register (~ A): 



A 



N>. 



In the more general case where the A^ are real num- 
bers with a finite resolvable phase resolution, a similar 
counting argument can be made as follows. Assume that 
the smallest measurable increment in phase Acj) is the 
same for all lattices, and divide the phase interval, to 
2tt, into bins of size A0. For the purpose of our scaling 
estimate (and unlike in the exact numerical calculation 
described above), phase is thus discretized, into 27r/A0 
bins which correspond to the approximate number of re- 
solvable states of each lattice ll|. The total number of 
states of all lattices is 



27r\ 
A^J 



The finest resolvable change in position is -Dmin = 
AiA(/)/27r, where Ai is the smallest lattice spacing. As the 
position is incremented in steps of -Dmin, the phases of all 
lattices necessarily go back to their original state within 
-^max steps. Hence the maximal representable range is 
bounded by 



Nx-1 



(3) 



As in standard RNS, the maximal combinatorial capacity 
may be achieved only for particular choices of A^. For 
generic choices of the lattice periods A^ , the representable 
range may not reach the bound of Eq. 0; nevertheless, 
we still expect the range -Dmax to scale exponentially with 
the number of lattices, and algebraically with the number 
of states in each lattice: 



Dmax OC ( 



(4) 



where the parameter a is expected to be of order unity. 
This relation is an ansatz, which we test by compari- 
son with the numerical results for particular sets of the 
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lattice spacings. We find that the maximum range typi- 
cally does scale according to the ansatz, as demonstrated 
Fig. 1 (solid lines in parts a and b). In part a, Dmax scales 
roughly exponentially with N\, and a ~ 0.55. In part b 
the dependence on phase resolution is shown for a partic- 
ular number of lattices, N\ = 12. In agreement with the 
ansatz, D-a^ay, scales roughly as a power of 27r/A0, with 
an exponent NeS c^i 9.7. 

In summary, our numerical results imply that capacity 
scales algebraically with phase resolution and exponen- 
tially with the number of lattices. Interestingly, this im- 
plies that to achieve a large capacity with a fixed number 
of neurons, neurons should be devoted to building more 



lattices rather than increasing the phase resolution: this 
result may help to explain the surprisingly poor phase 
resolution found in dMEC, where neural activity blobs 
cover 1 /3 of the total lattice period. 



[1] Note that the distance measure in Eq. Q, used for the 
exact numerical calculation, is written in terms of contin- 
uous phases and unlike in the scaling estimate here does 
not involve any binning. 



